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Abstract 


Theory indicates that pressure-temperature 
cross-spectra are necessary in predicting noise 
in regions of velocity gradients. In the present 
study, the feasibility of measuring pressure- 
temperature cross-spectra and coherence and 
temperature-temperature cross-spectra and coher- 
ence at spatially separated points along with 
pressure and temperature auto-spectra in a com- 
bustion rig was Investigated. The measurements 
were made near the inlet and exit of a 6.44 m long 
duct attached to a 0-47 combustor. For this study 
the fuel used was Jet A. The cross-spectra and 
coherence measurements show the pressure and tem- 
perature fluctuations correlate best at low fre- 
quencies. At the Inlet the phenomena controlling 
the phase relationship between pressure and tem- 
perature could not be Identified. However, at 
the duct exit the phase angle of the pressure is 
related to the phase angle of the temperature by 
the corrected flow time delay. 

Nomenclature 


Cq speed of sound 

c p specific heat at constant pressure 

f frequency 

G}-j thermocouple auto-spectrum 

G-jj thermocouple cross-spectrum 

H^ ratio of thermocouple auto-spectrum to 

cross-spectrum, G-j-j/G-n 
i vCT 

k wavenumber, w/c 0 

n integer 

p acoustic pressure 

R duct exit reflection factor 

S entropy perturbation 

T temperature perturbation 

W mass flow rate 

x displacement of point from combustor 

entrance 

y ratio of specific heats 

a square root of the variance of a 

distribution 

t thermocouple time constant 

u angular frequency, 2>rf 

Superscripts: 

* complex conjugate 

' perturbation 

Subscripts: 

g gas 

m measured 

1,2 thermocouple number 

3,4 duct locations 

Introduction 


A program to study turbofan engine core noise 
has been conducted at the NASA Lewis Research 


Center. In the study of core noise, the influence 
of combustion must be considered because the com- 
bustion process Is a major source of noise from a 
turbofan engine core. In a combustor, fluctua- 
tions in the heat release rate can cause both 
entropy ahd pressure perturbations. The pres- 
sure waves propagate through the engine and are 
radiated to the far field. The entropy variations 
are not direct noise sources because they have no 
associated pressure disturbance. However, as 
the entropy convects through regions of velocity 
gradients they, too, generate acoustic waves. 

The resulting combustion noise represents a 
combination of waves due to both sources. 

A model for predicting noise propagation in a 
variable area core nozzle which included convected 
entropy is discussed in Ref. 1. Using this model, 
calculations show that the noise level at a given 
frequency can be related to the entropy auto- 
spectrum. Furthermore, the core nozzle transfer 
function is affected by the entropy-pressure cross 
spectrum. The calculations show that if the 
entropy and pressure are in phase at the inlet, 
the magnitude and phase of the transfer function 
can have a rippled appearance. Consquently, in 
order to describe certain characteristics of com- 
bustion noise propagation, information about the 
cross-spectrum between the convected entropy and 
pressure and the entropy auto-spectrum are re- 
quired. However, no experimental entropy auto- 
spectra or entropy-pressure cross-spectra were 
available to use in the model. 

An entropy perturbation, S', is related to a 
perturbation in gas.pressure, p ' , and temperature, 
T' by 



(1 - r) P 1 
* Po 


( 1 ) 


Thus in order to obtain the entropy auto-spectrum 
and the entropy-pressure cross spectrum the tem- 
perature and pressure auto-spectra as well as 
the temperature-pressure cross-spectrum must be 
measured. The purpose of this paper is show that 
it is feasible to measure temperature-pressure 
cross-spectra in a combustion rig over a range 
of frequencies of interest in combustion noise 
studies, to obtain some baseline data, and to 
determine if some, non-random phase relationship, 
as hypothesized in Ref. 1, between pressure and 
temperature can exist. 

Dynamic temperature and pressure measurements 
have been made in combustors^ -5 in order to 
obtain auto-spectra. In this paper an exploratory 
study investigating the measurement inside a com- 
bustion rig of pressure-temperature cross-spectra 
and coherence at a single point as well as tem- 
perature cross-spectra and coherence at spatially 
separated points is presented. Pressure and tem- 
perature auto-spectra inside the combustion rig 
are also obtained. The engineering data acquired 


This paper Is declared a work of ihc U.S. 
Government and therefore Is in the public domain. 
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can be used In models, such as that discussed In 
Ref. 1, and compared with other published results. 


between the measured and theoretical values of 
Hjj and H22 where 


The measurements to be described were made in 
a long (6.44 m) constant area duct attached to a 
0-47 combustor. The rig was designed to operate 
with no combustion in the long duet. However, for 
all the test conditions some combustion did occur 
at the duct inlet station where one set of pres- 
sure and temperature measurements were made. 
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1 + W T 


TT 


and 


( 2 ) 


Experimental Apparatus 

The temperature and pressure measurements were 
made in the combustion r'ig shown in Fig. 1 which 
consisted of a J— 47 can combustor (see Fig. 1(b)) 
mounted in a 0.3 m diameter by 0.77 m long test 
section and an attached 6.44 m long 0.3 m diameter 
stainless steel duct. The fuel used was Oet A. 
This paper discusses pressure and temperature mea- 
surments at stations 3 and 4. The test conditions 
are shown in Table I. 

All tests were conducted in an outdoor acous- 
tic arena. The facility is shown in Fig. 2. 

Note that in all cases the duct exit tempera- 
ture was greater than the duct inlet temperature 
indicating some combustion was occuring in the 
duct. This was due to mechanical problems which 
developed iri the combustor during the test pro- 
gram. These problems make the test conditions 
difficult to duplicate and hinder the inter- 
pretation of the measurements. However, the 
combustor mechanical problems have no bearing on 
the examination of the feasibility of measuring 
temperature-pressure cross spectra. 

Pressure Measurements 
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where 


i * V-I 


The temperature auto-spectra are Gjj, and G22 
and the temperature-temperature cross spectrum is 
G12 • The model used assumes that the measured 
temperature, T m , is related to the gas tempera- 
ture, T g by 


1 + 1 ® 


(4) 


in the frequency domain. Futhermore, it is assumed 
that the coherence between the two thermocouples 
is high so that extraneous noise can be neglected. 

Examination of the plots of the real and 
imaginary values of H-j-j can yield information on 
the time constants since the imaginary value of 
Hi i has a maximum at id = 1/x-j. Furthermore 


The internal pressure transducers used were 
conventional 0.635 cm microphones with pressure 
response cartridges. In order to avoid direct 
exposure of the microphone to the combustion 
gases, they were mounted outside the ducted com- 
bustion rig and the fluctuating pressure in the 
rig was transmitted to the transducers by means of 
a ’‘semi-infinite" acoustic wave guide. Details on 
these probes are given in Ref. 6. 

Temperature Measurements 


11m H^ = tg/Tj. ( 5 ) 

(jJ+CO 

and 

1 1m ^22 c t i / t 2 (6) 

(0-foo 

The time constant in all cases was determined 
using the spectrum between 0 and 40 Hz. 


The temperature fluctuations were measured 
using chromel-alumel thermocouples. The two 
wire thermocouple probe design shown in Fig. 3 
was selected so’that the response time of the 
thermocouple could be determined based on the 
actual test data, as opposed to a separate pre- 
calibration. The procedure discussed in Ref. 7 
was used. The thermocouple chromel-alumel wire 
diameters were 76.2 and 25.4 pm. The thermo- 
couples were inserted into the duct by actuators 
after combustion started in order to prevent their 
destruction by the combustion start up transient. 
While the thermocouples were quite durable, their 
lifetime was drastically reduced by operation at 
elevated temperatures (1115 K). 

The approach used to determine the themocouple 
time constants is based on the two thermocouple 
method described in Ref. 7 However, the actual 
time constants were not found by examination of 
the plots as suggested in Ref. 7. Instead, a non- 
linear minimization technique was used to find the 
time constants to minimize the least square error 


In order to demonstrate the procedure, typical 
constant bandwidth (0.08 Hz) two wire thermocouple 
auto-spectra, cross-spectra and coherence function 
measurements are shown in Figs. 4 to 8. Note that 
the thermocouple coherence is high (Fig. 6). This 
is a necessary condition for application of the 
calibration procedure since it means that extrane- 
ous noise is low. In addition, the measured 
values of Hji and H22 are shown in terms of 
their real ana imaginary components along with the 
least square curve fits used to find the time con- 
tants (Figs. 7 and 8). The <u = 1/x-j points and 
the high frequency asymptotic value of H-j-j are 
also identified. Also, in Fig. 4 the gas tempera- 
ture, T g , auto-spectrum and cross-spectrum after 
correction are shown. Table I gives the time con- 
stants used to correct the cases used herein. 

Results and Discussion 

This section discusses measurements of 
pressure-pressure and temperature-temperature 
cross-spectra between the inlet and the exit of 
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the long duct. Also, pressure-temperature cross- 
spectra at the Inlet and the exit of the duct are 
discussed. Only, a few typical results sufficient 
to establish definite trends are Included here. 

With each cross-spectrum, the corresponding 
coherence function is presented. Note that a 
coherence smaller than unity implies one or more 
of the following condi tians exists: 

1. Extraneous noise is present in the mea- 
surements. 

2. The system relating the output to the input 
Is non-linear. 

3. The output is due to more than one Input. 
Pressure Cross-Spectra and Auto-Spectra 

Typical constant bandwidth (2 Hz) pressure to 
pressure cross-spectra and pressure auto-spectra 
are shown In Figs. 9 to 11. The cross-spectra 
magnitude plots and auto-spectra magnitude plots 
have certain slmllarltes since as the temperature 
increases, the frequency at which a particular 
spectral feature occurs increases in a proportion 
that keeps the wave number constant. However, the 
variation in temperature along the duct produces 
distortions In the shape of the magnitude plot 
that would not appear in a constant temperature 
duct. The pressure auto-spectra are shown so that 
they can be compared with the corresponding 
temperature-pressure cross-spectra. 

The cross-spectra phase relationships are due 
to a longitudinal standing wave in the duct. 

The cross-spectra phase angle Is inti al ly zero 
corresponding to a positive value of P3P4. The 
phase angle switches by 180 degrees at approxi- 
mately 50 Hz intervals. The phase angle variaton 
with frequency of the pressure cross spectra In 
the duct Is understandable In terms of a simple no 
flow model for the acoustic pressure dlstrlbtlon 
In the duct. Without flow the pressure in the 
duct is proportional to 


Furthermore, at low frequencies the reflection 
factor is nearly - 1, Consequently, the pressure 
is approximately proportional to sin(kx ). Thus 
the pressure cross-spectrum between station 3 and 
4, P3P4* in a duct without flow is then propor- 
tional to 

p 3 p* * s1n(kx 3 )sin{kx 4 ) (8) 

where X3 = 1.47 m and X4 ■» 6.91 m. This func- 
tion Is either positive or negative and it changes 
sign at frequencies n.f which are multiples of 
the first zero of P3P4. The phase angle devia- 
tions from this simple model are due to the pres- 
ence of flow and the duct reflection and the duct 
reflection factor not being -1. 

The corresponding coherence functions are 
shown in Fig. 12. The coherences are quite large 
above 20 Hz. However, at multiples of the first 
zero of P3P4 the coherence is zero. 

Temperature-Temperature Cross-Spectra 

Temperature-temperature constant bandwidth 
(0.4 Hz) cross-spectra between the entrance and 
the exit of the duct are shown for three cases in 


Fig. 13. All temperature measurements are shown 
after correction The cross-spectrum magnitude is 
largest at low frequencies. 

At low frequencies (up to 50 Hz) the phase of 
the cross spectra changes linearly with frequency 
with slopes proportional to the time delay due to 
convection of temperature In the duct. The time 
delay 1$ about 125 ms for the data shown In Figs. 
13(a) and (b) where the mass flow rate, W, ,is 
0.5 kg/s and 62.5 ms for the data shown In Fig. 
13(c) where the mass flow rate Is 1.13 kg/s, 

The corresponding coherence data are shown in 
Fig. 14. At the higher frequencies the coherence 
is very low (Fig. 14) and no orderly relationship 
between the temperatures is observed. Note also, 
from Fig. 14, that as the amount of burning In the 
duct increases the coherence is greatly reduced. 

The thermocouple signals show the presence of 
60 Hz and 110 Hz noise. The signal at station 4 
at the duct exit was more contaminated than the 
signal at station 3 at the duct inlet, 

Pressure-Temperature Cross-Spectra 

Duct inlet : Pressure-temperature constant 
bandwidth (0.40 Hz) cross-spectra at the duct 
inlet (station 3) for three cases are shown in 
Fig. 15. The corresponding coherence fuetlons 
are shown in Fig. 16. The mechanical cut-off fre- 
quency of the microphone due to pressure equaliza- 
tion is about 5 Hz. Consequently, below this 
frequency the data are unreliable. 

Above 50 Hz the shapes of the pressure- 
temperature cross-spectra magnitude plots (Fig. 

15) resembles the shapes of the corresponding 
pressure auto-spectra magnitude plots shown In 
Fig. 10. This Indicates that the part of the 
temperature perturbation coherent with the pres- 
sure perturbation has no strong spectral char- 
acteristics above 50 Hz. 

The phase angle variation shows distinct 
atterns, indicating a systematic relationship 
etween pressure ana temperature near the source 
region. However, the physical phenomena producing 
this pattern have not been identified. 

Comparing, the pressure-temperature cross- 
spectrum magnitude plots (Fig. 15) and coherence 
plots (Fie* 16 ) with the pressure auto-spectra 
plots (Fig. 10) shows that the temperature per- 
turbation is weakly excited in certain frequency 
bands near particular duct resonant frequencies. 
For inlet temperatures of 1115, 875 and 920 K note 
the 110 Hz, 98 and 70 Hz regions respectively. 

The systematic relationship between the pres- 
sure and the temperature at the inlet is not a 
strong one. Consequently, it is more visible in 
the phase plots than the coherence plots. 

Duct exit : Pressure-temperature constant 
bandwidth (0.40 Hz) cross-spectra at the duct 
exit (station 4) for three cases are shown in 
Fig. 17, The corresponding coherence fuetlons are 
shown in Fig. 18. Again the magnitude plots of 
the pressure-temperature cross spectra resemble 
the magnitude plots of the corresponding pressure 
auto-spectra shown in Fig. 11. 
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The phase of the cross spectra again show 
linear dependence with slopes that are propor- 
tional to the time delay due to convection of 
temperature In the long duct. Note, that since 
this is a phase measurement between two different 
variables, that Is pressure and temperature, at 
the same point, the time delay must be due to the 
difference In travel times between a pressure and 
temperature disturbance as they travel from the 
same upstream source region. This shows that tem- 
perature and pressure fluctuations are related at 
the combustion source point. Furthermore, an 
extrapolation of the cross-spectrum suggests that 
the two are in phase at low frequencies. 

In view of these findings that pressure and 
temperature are in In phase at the source in this 
combustion rig, it is possible that the entropy 
and pressure in the JT15D nozzle might also be in 
phase at the nozzle inlet. Thus, some of the 
ripple in the measured JT15D transfer function 
presented in Ref. 1 might be due to convected 
entropy. This suggests that it might be of 
interest to obtain dynamic temperature measure- 
ments along with acoustic pressure measurements in 
turbofan engine tests. 

Temperature Auto-Spectra 

Because plots of the gas temperature auto- 
spectra resemble those presented in Fig. 5, addi- 
tional ones are not included. It can be seen that 
most of the fluctuations are at low frequencies. 
The temperature standard deviations calculated by 
integrating the temperature auto-spectra for a 
frequency range from 0 to 200 Hz are shown in 
Table I. The percent r.m.s. temperature fluctua- 
tions are also shown in Table I. At the inlet 
where combustion is occuring they vary in a range 
from 15 to 20 percent. At the duct exit the range 
is 2 to 4 percent. These results are similar to 
those found at the exit of the gas turbine combus- 
for by Dils 1 2 where the per cent r.m.s. temperature 
fluctuations were 7 and 10 percent at a takeoff 
operating condition where the average temperatures 
were 1430 and 1438 K. The burner exit r.m.s. 
temperature fluctuations reported in Ref. 5 are 
generally lower than those observed by Dils. 

Conclusions 

From the results of this exploratory study, 
the following conclusions were reached: 

1. Pressure and 'temperature cross-spectra and 

coherence functions can be measured in a combus- 
tion rig using a two-wire thermocouple probe and 
an infinite tube microphone over a range of fre- 
quencies of interest in combustion noise studies. 


2. The temperature and pressure fluctuations 
are related at the combustion source point and an 
extrapolation of the cross-spectrum phase suggests 
that the two are in phase at low frequencies. 

3. At the duct exit the phase angle of the 
pressure is related to the phase angle of the 
temperature by the convected flow time delay. 

These three conclusions, taken together, 
provide confidence that identifiable physical 
relationships between temperature and pressure 
in a combustion system can be measured using 
sophisticated data analysis techniques. These 
relationships are valuable either as inputs or 
as verification of outputs of theoretical 
models describina combustion noise generation 
or propagation. 
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TABLE I, - THERMOCOUPLE TIME CONSTANTS, TEMPERATURE 
STANDARD DEVIATION, AND RMS TEMPERATURE FLUCTUATION 


Test Point 

87 

86 

88 

T exit- K 

1135 

990 

1120 

K 

1115 

875 

920 

W, kg/s 
Time constant 
Station 4 (exit) 

0.5 

0.5 

1.13 

(76.2 pm dlam), ms 

41 

79 

49 

(25.4 pm dlam), ms 
Station 3 (Inlet) 

20 

34 

22 

(76.2 pm diam), ms 

61 

67 

42 

t 2 (25.4 pm dlam), ms 
Standard temperature deviation, o 

36 

40 

25 

Station 4 (exit), K 

35 

40 

25 

Station 3 (inlet), K 
Percent rms temperature fluctuation 

170 

150 

170 

Station 4 (exit) 

3 

4 

2 

Station 3 (inlet) 

15 

17 

19 



X TWO WIRE PfiOBE 
(a) Rig schematic. 

Figure 1. - Schematic of ducted combustion system rig. 












Figure 2. - NASA Lewis combustion acoustics facility. 
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Figure 4 - Two wire chromel alumel thermocouple time 
constant determination; temperature cross-spectrum 
between thick and thin thermocouples, G 12 (station 4, 
Tmijrr = 875 K, Tgxfx ” 990 K, W “ 0. 5 kg/s, band- 
width - 0. 08 Hz). 
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(a) 76. 2 |im diameter (G^. 

(b) 25, 4 pm diameter (G^). 

Figure 5. - Two wire chromel alumel thermocouple 
time constant determination: temperature auto- 
spectrum of thick and thin thermocouples (station 
4, Tiwirr c 875 K, W “ 990 K, W ■ 0.5 kg/s, 
bandwidth - 0.08 Hz). 



Figure 6. - Two wire chrome! alumel thermocouple 
time constant determination*, temperature coherence 
between thick and thin thermocouples (station 4, 
t INLET e 875 K . Texit - 990 K , W - 0. 5 kg/s . band- 
width - 0.08 Hz). 
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FREQUENCY, Hz 

(a) T|nue T " 1H5 K, T EX | T " 1135 K, W - 0.5 kg/s. 

Figure 9. - Pressure cross-spectra between duct inlet 
(station 3} and duct exit (station 4), P 4 P 3 * (bandwidth 
* 2 Hz). 
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FREQUENCY, Hz 


W ^INLET " 875 K, T EX |j - 990 K, W - 0.5 kg/s. 
(c ^ t INL£T“ 920 ’EXIT “ H20 K, W- 1.13 kg/s. 


Figure 9. - Concluded. 
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FREQUENCY, Hz 

(a) T | NLEr «m5K. T ex , t - 1135 K, W- 0.5 kg/s. 

Figure 13. - Temperature cross-spectra between duct inlet 
(station 3) and duct exit (station 4), T 4 T 3 * (bandwidth ■ 
0.4 Hz). 











FREQUENCY, Hz 

^INLET B K* t EXU b H35 K, W ■= 0.5 kg/s. 

Figure 14. - Temperature coherence between duct inlet (sta- 
tion 3) and duct exit (station 4), I TJo* I / i T,, I 2 1 To |2 
(bandwidth ■ 0. 4 Hz). y J 





0 32 64 96 128 160 


FREQUENCY, Hz j 

(b) Tj NLET - 875 K, T ex , t - 990 K, W ■ 0.5 kg/s. 1 

(c) T INLET * 920 K, T EX j T - 1120 K, W» 1.13 kg/s. | 

Figure 14. - Concluded. j 



PHASE deg MAGNITUDE, dB, (Po«20tiPa 

T 0 =44jiK) 
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Figure 15. - Cross-spectra between pressure and temperature 
at duct Inlet (station 3), p 3 T 3 * (bandwidth - 0.4 Hz). 
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0 32 64 96 1 28 160 


FREQUENCY, Hz 

(a,T INLET = 1115K ' t EXIT " 1135 K * W‘ 0.5 kg/s. 

(b) T|^LET c **75 K, Texit B 990 K, W c 0.5 kg/s. 

(c) t INLET= 920 K, T EX | T »1120K, W = 1.13 kg/s. 

Figure 16. - Coherence between pressure and temperature 
at duct inlet (station 3), I p 3 T 3 * I 2 / 1 T 3 I 2 1 p 3 1 2 (band- 
width b 0. 4 Hz). 
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FREQUENCY, Hz 


(a) T wlET . ms k, T EX | T - 1135 K, W - 0. 5 kg/s. 

Figure 17. - Cross-spectra between pressure and temperature 
at duct exit (station 4), (bandwidth - 0.4 Hz). 
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Figure 18. - Coherence between pressure and temperature 
at duct exit (station 4), i p 4 T/ I 2 / IT 4 I 2 I p 4 1 2 (band- 
width - 0.4 Hz). 



FREQUENCY, Hz j 

<&)T| N LET - 875 K, T EX|T - 990 K, W - 0.5 kg/s. * j 

(0 T| NLEr ■ 920 K, T EX)T - 1120 K, W - 1. 13 kg/s. j 

Figure 18. - Concluded. ’ j 
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